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PREFACE 

This  Memorandum  was  prepared  as  part  of  Pj\ND's  continuing 
study  of  advanced  radar  signal,  processing  methods.  It  extends 
previous  work  on  optimum  signal  processing  methods  for  determining 
the  characteristics  of  a  radar  target  by  analyzing  the  target's 
radar  signature.  The  results  of  the  study  are  applicable  to  the 
analysis  of  earth  satellites  and,  suitably  extended,  could  be 
applied  to  radar  astronomy  and  radar  reconnaissance. 
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SUMMARY 


A  problem  of  considerable  interest  in  a  variety  of  applications 
is  that  of  estimating  the  characteristics  of  a  radar  target  by 
analyzing  its  radar  signature.  Applications  of  interest  include 
satellite  identification,  radar  astronomy,  and  radar  ground  mapping. 
Target  characteristics  of  interest  include  not  just  range,  range  rate, 
and  range  acceleration,  but  also  dimensions,  shape,  and  rate  and  axis 
of  rotation. 

Although  many  techniques  have  been  developed  to  attack  these 
problems,  the  theory  of  optimum  signal  processing  methods  and  of  the 
optimum  estimation  accuracies  to  be  expected,  as  a  function  of  the 
noise  level,  the  data  set,  and  the  a  priori  knowledge  of  the  target 
parameters,  is  still  not  well  developed. 

This  Memorandum  extends  the  inventory  of  techniques  for  defining 
optimum  signal  processing  methods  and  predicting  the  accuracies 
achieved.  Chief  attention  is  given  to  estimation  of  linear  or  quad¬ 
ratic  functionals  of  the  scatterer  amplitudes  (such  as  the  average 
scattering  cross  section  attributable  to  specific  regions  on  the 
target);  however,  the  estimation  of  other  target  parameters  is  also 
treated.  The  analysis  is  carried  farthest  for  the  signature  analysis 


of  rotating  targets. 
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I .  INTRODUCTION 

The  analysis  reported  in  this  Memorandum  was  initiated  to 
develop  a  theory  which  would  predict  the  accuracy  with  which  the 
characteristics  of  an  extended  radar  target  can  be  determined  by 
analyzing  its  radar  signature.  "Extended"  means  extended  compared 
to  a  wavelength,  not  necessarily  compared  to  a  radar  resolution  cell. 

In  fact,  some  of  the  interesting  cases  are  those  in  which  the  target 
dimensions  are  much  smaller  than  a  spatial  resolution  cell  of  the 
radar.  Of  primary  interest  are  rotating  targets,  although  considera¬ 
tion  will  not  be  entirely  restricted  to  this  case.  Characteristics 
of  interest  to  determine  might  include  size,  shape,  rate  of  rotation, 
axis  of  rotation,  or  radar  cross  section  of  specified  portions  of  the 
body.  This  problem  is  of  interest  in  identification  of  earth 
satellites  and  in  radar  astronomy,  where  the  target  is  extended  but 
isolated;  and  also  in  ground  mapping,  where  the  target  is  an  extended 
continuum.  (In  ground  mapping,  and  in  signature  analysis  of  stabilized 
satellites,  the  targets  are  not,  strictly  speaking,  rotating.  However, 
even  in  these  cases,  the  relative  motion  of  the  target  with  respect 
to  the  radar  can  often  be  considered,  to  a  good  approximation,  a 
rotation  about  the  line  of  sight,  after  mean  translatory  motion  has 
been  corrected  for.)  Another  type  of  problem  which  can  be  considered 
within  this  framework  is  the  estimation,  in  the  presence  of  noise,  of 
the  illumination  pattern  of  a  rotating  antenna  by  observations  of  its 
far  field. 

A  great  deal  of  literature  exists  on  this  subject.  Many  special 
methods  have  been  evolved,  and  a  certain  amount  of  theory  has  been 
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developed  to  predict  the  effects  of  various  relevant  factors  on  the 
success  which  can  be  obtained.  Nevertheless,  it  cannot  be  said  that 
a  satisfactory  theory  has  been  developed.  By  a  satisfactory  theory, 
is  meant  one  which  will  specify  good  or  optimum  signal  processing 
methods  to  determine  target  characteristics,  and  will  also  predict 
the  accuracy  or  reliability  with  which  target  characteristics  of 
interest  can  be  estimated,  as  a  function  of 

a.  types  and  statistical  characteristics  of  errors, 

b.  characteristics  of  the  data  set,  and 

c.  degree  of  a  priori  knowledge  of  the  target. 

Types  of  errors  might  include,  for  example,  additive  noise 
errors,  phase  errors,  multiplicative  amplitude  noise,  nonlinear 
distortions,  etc.  Characteristics  of  the  data  set  would  include 
waveform  design,  types(s)  of  radar  (monostatic,  multistatic,  inter¬ 
ferometric),  extent  of  frequency  ar J  polarization  diversity, 
scattering  matrix  elements  to  be  measured,  degree  of  coherence  (not 
only  temporally,  at  one  sensor  location,  but  from  sensor  to  sensor), 
and  length  of  observation  interval  or  intervals.  A  priori  knowledge 
about  the  target  might  take  the  form,  to  cite  just  a  few  examples, 
of  knowledge  that  it  is  (from  the  electromagnetic  viewpoint)  a 
spherical  surface,  as  in  radar  astronomy;  that  it  is  a  stabilized 
body;  or  that  it  is  a  long,  slender  body. 

At  present,  there  are  many  relatively  simple  questions,  in 
relatively  simple  contexts,  which  cannot  be  satisfactorily  answered 
by  existing  theory.  An  example  closely  related  to  the  situations 
which  this  Memorandum  analyzes  is  the  "low  resolution"  signature 
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analysis  of  radar  targets  such  as  earth  satellites.  It  is 
well  known  that  various  characteristics  of  a  target  of  this  type, 
such  as  dimensions  and  to  some  extent  shape,  can  be  estimated  rather 
well  if  the  target  can  be  viewed  over  a  reasonable  set  of  angular 
aspects,  preferably  but  not  necessarily  including  the  broadside  aspect; 
and  that  such  estimation  can  be  performed  even  if  the  spatial  radar 
resolution  cell  is  much  larger  than  the  physical  dimensions  of  the 
target.  Moreover,  this  does  not  necessarily  require  phase-coherent 
methods.  However,  there  is  at  present  no  theory  which  will  predict 
how  the  ability  to  estimate  specified  target  characteristics  depends 
on  the  noise  level  and  on  the  set  of  aspects  over  which  the  body  is 
observed.  The  effect  of  these  factors  must  be  determined  empirically, 
by  submitting  signature  records  to  analysts. 

A  satisfactory  theory  must  also  have  the  property  that  answers 
may  be  obtained  in  reasonable  time  and  at  reasonable  cost  on  digital 
computers . 

This  Memorandum  by  no  means  succeeds  in  establishing  a  fully 
satisfactory  general  theory  of  target  identification  by  radar  signature 
analysis.  In  fact,  a  look  at  the  problem  rather  quickly  reveals  that 
it  is  difficult  to  identify  problems  which  are  based  on  reasonable 
models  of  situations  of  practical  interest  and  which,  at  the  same  time, 
can  be  solved,  even  using  digital  computers.  The  objective  of  this 
study  is  to  contribute  to  the  methodology  available  for  attacking 
these  problems  by  considering  examples  based  on  models  reasonably 
near  to  physical  reality--at  least  in  the  sense  that  the  results 
obtained  would  give  considerable  insight  into  realistic  situations-- 
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and  for  which  it  is  feasible  to  carry  out  the  required  calculations. 
Thus,  what  is  derived  here  is  not  a  general  theory,  but  it  is  an 
addition  to  the  arsenal  of  special  theories  for  this  class  of  problems. 
The  actual  implementation  of  the  indicated  calculations  requires 
(except  in  very  special  cases)  computer  programs  for  solving  certain 
integral  equations;  this  has  not  been  done  as  yet,  although  it  will 
be  clear  that  the  necessary  computations  are  well  within  the  realm  of 
feasibility . 

POSSIBLE  APPROACHES 

Theoretical  approaches  to  the  class  of  problems  just  stated 
include  the  following: 

o  Modeling  the  electromagnetic  signal  scattered  from  the 
extended  target  as  a  noise  process,  and  estimating  the 
parameters  of  the  noise  statistics  (swch  as  parameters  of 
the  covariance  function).  Reference  1  provides  good  examples 
and  a  reasonably  comprehensive  bibliography.  Essentially, 
this  amounts  to  a  model  of  the  target  as  a  reasonably  numerous 
collection  of  independent  scatterers. 

o  Modeling  the  scattered  signal  as  depending  predominantly  on 

a  finite  set  of  parameters  which  can  be  identified  with 

characteristics  of  the  target  motion  or  configuration,  and 

applying  parameter  estimation  theory,  usually  based  on 

maximum-likelihood  methods  and  Cramer-Rao  bounds  (equivalently, 

"information  matrices")  for  estimating  attainable 

(2  3) 

accuracies.  ’  For  example,  the  target  might  be  considered 
to  consist  of  some  reasonably  small  number  of  rigidly 
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connected  point  scatterers;  the  unknown  parameters  would 
then  include  the  radar  cross  sections  and  relative  positions 
of  these  scatte  ers . 

Although  many  useful  results  have  been  obtained  by  these  approaches, 
both  have  drawbacks  when  applied  to  the  problem  at  hand.  While  the 
treatment  of  the  signal  as  a  noise  process  does  not  absolutely  require 
the  assumption  that  the  received  signal  process  is  Gaussian  or 
stationary,  those  treatments  which  attempt  to  define  optimum  signal 
processing  methods,  and  to  determine  lower  bounds  on  the  estimation 
errors,  usually  assume  that  the  noise  is  Gaussian,  and  often  assume 
stationarity  or  quasi-stationarity  in  some  sense. 

However,  many  examples  of  the  class  of  problems  under  considera¬ 
tion  are  characterized  by  the  fact  that  the  signal  statistics  are  not 
Gaussian,  or,  even  if  they  can  reasonably  be  modeled  as  Gaussian,  they 
are  highly  nonstationary.  Moreover,  this  non-Gaussian  or  non¬ 
stationary  behavior  often  involves  essential  aspects  of  the  problem. 

For  example,  the  signal  may  have  a  significant  or  dominating 
component  arising  from  a  very  small  number  of  scatterers.  Or,  more 
important  for  the  examples  to  which  this  Memorandum  is  mainly 
devoted,  the  scatterer  may  act  essentially  as  a  line  source  or  a 
plane  source,  i.e.,  there  may  be  a  small  set  of  aspects  (perhaps  only 
one  or  two)  at  which  all  scatterers  are  in  phase,  producing  what  is 
called  a  "specular  flash";  this  specular  flash  is  of  great  importance 
in  the  problem  of  estimating  target  characteristics. 

The  parameter  estimation  approach,  insofar  as  it  is  based  on 
Cramer-Rao  bounds,  information  matrices,  and  maximum  likelihood 
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methods,  has  the  drawback  that  it  is  often  difficult,  in  cases  of 
nonlinear  dependence  on  the  parameters  to  be  estimated,  to  say  whether 
the  accuracy  bounds  thus  obtained  are  informative;  that  is,  it  is 
difficult  to  determine  whether  the  Cramer-Rao  lower  bounds  on  error 
variance  are  reasonably  close  to  the  actual  minimum  attainable 
variances.  In  regular  estimation  cases,  the  Cramer-Rao  variance  bounds 
are  asymptotically  the  minimum  variances;  "asymptotically"  means  for 
sufficiently  large  signal -to-noise  ratio  (SNR).  However,  as  pointed 
out  in  Ref.  3,  there  are  many  pitfalls  in  determining  how  large  the 
SNR  must  be  for  the  Cramer-Rao  bounds  to  be  reasonably  close  to  the 
minimum  variances  even  if  the  estimation  problem  is  "regular."  It  is 
definitely  not  true  that  output  SNR  which  is  sufficient  for  reliable 
target  detection  is  sufficiently  large  to  be  in  the  asymptotic  region 
just  referred  to.  In  fact,  there  are  relatively  simple,  practical, 
and  regular  parameter  estimation  problems  (a  small,  highly  non- 
exhaustive  set  of  examples  is  given  in  Ref.  3)  in  which  the  Cramer-Rao 
bound  may  be  highly  uninformative  for  SNRs  which  are  adequate  for 
detection  purposes  (i.e.,  the  bound  may  be  a  small  fraction  of  minimum 
attainable  variance) .  This  difficulty  is  even  worse  in  more  complex 
estimation  problems,  such  as  those  involving  the  estimation  of  con¬ 
figuration  parameters  of  connected  sets  of  point  scatterers. 

The  Bayes  or  Barankin  approaches  to  parameter  estimation  avoid 
these  pitfalls  in  theory,  but  ordinarily  the  computational  difficulties 
associated  with  them  are  formidable.  Nevertheless,  the  approach 
selected  here  is  a  Bayes  approach.  This  approach  is  based  on  the 
minimum  mean  square  error  criterion.  Theoretically,  results  can  be 
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obtained  for  any  signal- to-noise  ratio.  This  is  also  practical  in  a 
number  of  cases  of  interest.  In  the  interest  of  computational 
feasibility)  hovirever,  it  will  in  other  cases  be  necessary  to  reintroduce 
some  form  of  strong-signal  assumption.  The  Bayes  approach  also  has 
the  advantage,  over  that  of  maximum  likelihood,  that  it  is  possible 
to  deal  directly  with  functions  of  the  unknown  parameters;  usually, 
functions  of  the  parameters,  rather  than  individual  parameters,  are 
of  greatest  interest. 

Since  a  priori  statistics  must  be  associated  with  the  unknown 
parameters  in  a  Bayes  approach,  it  is  also  possible  to  say  that  the 
scattered  signal  is  considered  to  be  a  random  process.  Although  most 
of  the  results  will  be  derived  under  assumptions  leading  co  a  priori 
Gaussian  statistics  of  the  signal,  some  of  the  results  can  be  extended 
to  cases  where  the  a  priori  signal  statistics  are  mixtures  of  Gaussian 
statistics,  i.e.,  they  result  from  choosing  a  Gaussian  probability 
measure  from  a  specified  set  of  such  measures,  with  a  probability 
distribution  over  the  set  of  possible  Gaussian  measures.  Also,  it  is 
possible  to  deal  with  cases  in  which  there  are  peculiar  forms  of  non- 
stationarity,  e.g.,  those  in  which  specular  flashes  occur.  In  addition, 
the  models  to  be  analyzed  also  are  capable  of  reflecting  various  degrees 
of  a  priori  knowledge  of  target  characteristics. 

Section  II  contains  the  derivation  of  the  general  methods  to  be 
employed.  The  target  models  under  consideration  in  that  section  are 
not  confined  to  rotating  targets.  Section  III  states,  for  a  selection 
of  rotating  target  models,  the  form  taken  by  the  more  general  expressions 


of  Section  II . 
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II.  GENERAL  RESULTS 

A.  INITIAL  FORMULATION  AND  RESULTS 

In  the  present  section,  the  class  of  problems  to  be  considered 
will  not  be  restricted  to  rotating  cargets.  Initially,  it  is  supposed 
that  the  received  signal,  in  the  presence  of  additive  noise,  can  be 
represented  by 

S(t)  =  ^a(x(i))  F (t,  x(i))  +  e(t)  (1) 

i 

where  e(t)  is  additive  noise.  S(t)  is  observed  over  a  known  time 
interval  or  intervals. 

Here,  a  is  a  function  of  the  argument  x,  which  in  turn  is  a 
possibly  multicomponent  vector,  ranging  over  a  finite  set  of  known 
vectors  {x^^},  i  =  1,  ...,  n.  This  could  be  regarded  as  an  approxi¬ 
mation  to  the  representation  of  S(t)  as  an  integral,  although  such  an 
interpretation  is  not  always  necessary.  If  such  an  interpretation  is 
adopted,  the  required  limiting  procedures  must  be  handled  with  great 
care  (see  Section  II. C). 

The  parameter  t  represents  time,  although  here,  too,  it  is 
possible  to  consider  t  to  be  a  multicomponent  vector  (although  this 
will  not  be  reflected  in  the  notation).  For  example,  if  there  are 
K  different  sensors,  it  would  be  possible  to  have  t  =  (t^,  ..., 
with  t  ranging  over  the  observation  times  of  the  k*"*1  sensor. 

K 

Similarly,  if  a  given  sensor  observes  several  different  cross¬ 
polarizations,  this  could  be  formulated  by  considering  t  to  be  multi¬ 
dimensional.  The  vector  £  will  represent  primarily  position  co¬ 
ordinates  with  respect  to  a  coordinate  system  fixed  to  the  target. 
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Other  components  may  be  adjoined  to  x  to  represent  internal  velocity 
components  of  the  target,  if  any,  and  also  to  differentiate  between 
two  phase  components  of  the  signal  scattered  from  any  given  point  on 
the  body.  This  will  all  be  illustrated  by  concrete  examples  in  the 
next  section. 

The  set  of  unknown  parameters  is  a(x^),  i  =  1,  2,  ...,  n. 

The  problem  which  will  primarily  be  addressed  is  the  estimation  of 
real-valued  functionals  of  a,  and  most  important,  of  linear  and 
quadratic  functionals  of  a. 

It  is  assumed  that  the  a  priori  statistics  of  both  {aQO}  and 

£  e  ( t ) }  are  Gaussian,  and  that  (cKx)}  is  statistically  independent  of 

[e(t)}.  The  assumption  of  Gaussian  a  priori  statistics  makes 

S(t,  oi)  a  Gaussian  process;  later,  cases  will  be  considered  in  which 

F  depends  on  additional  unknown  parameters  X,  in  which  cases  S  would 

* 

no  longer  be  a  Gaussian  process. 

It  will  further  be  assumed,  for  expository  purposes,  that 
[a(x)}  and  [e(t)}  have  zero  ensemble  means,  and  that 


E 


<y(x(i))  a(x(j)) 
E  j^e(t)  e(t') 

E  |_a(£(i))  e(t) 


A.  . 
ij 


<P(t,  t') 

0 


(2) 

(3) 

(4) 


A  .  . 

The  assumptions  about  the  statistics  of  £cy  (x) }  are  assumptions 
about  the  a  priori  state  of  knowledge  used  in  forming  Bayes  estimates, 
and  do  not  necessarily  apply  to  any  particular  scattering  target. 
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It  is  then  possible  to  state  the  conditional  (a  posteriori) 
probability  density  function  (p.d.f.)  of  a  given  S: 

_  n  \ 

p(a  I  §,)  =  (2n)2  j^det  B~j 

x  exe  {-*  I  Bij  [“<£(i)>  -  S(5(i))]  [a(x(j))  -  5(x(Jh] 


J 


where 


Q'(x(l))  =  lim  £  Y  B-j  F(t  ,  x(j))  S (t  ) 


M-jU  j 

5.  .  =  A?1  +  lim  7  T]  F(t  ,  xv"')  F(r  xVJ') 
ij  ij  L  Vu  \i  -  u  ~ 


(i>N  r./.  „(j)' 


th 


T]  r  =  (p,,  u)  element  of  matrix  inverse  to  qj ( t  ,  t^) 


‘p,u 


u  u' 


In  these  formulas,  lim  implies  the  limit  as  the  points  tt^) 
become  dense  in  the  observation  interval.  Such  limits  are  rigorously 
defined  if  (e(t)}  is  continuous  in  the  mean  and  the  functions 

/  •  \  / o\ 

F(t,  x/  ;)  satisfy  certain  regularity  conditions.  In  Eqs .  (6)  and 

(7),  B_  and  refer  to  the  (i,  j)  elements  of  the  matrix  inverses 
of  B  and  A.  There  is  a  large  body  of  literature  on  methods  of 

evaluating  expressions  of  the  type  appearing  in  Eqs.  (6)  and  (7); 

Ref.  3  contains  a  partial  summary  of  such  methods. 


It  will  prove  convenient  to  write  Eqs. 
equations  below,  in  vector-matrix  notation: 


(5) 


(6) 


(7) 


(8) 


a  =  B_1  FT  T)  S 

***  v«.  *** 


-1  T 

B  =  A  +  F  T)  F 


(9) 

(10) 


(6)  and  (7),  and  similar 
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a  ^  ( i ) 

Here,  a  is  an  n-component  vector  with  elements  o' (x  );  A  and  B 

MM  *** 

are  n  x  n  matrices;  and 


(FT  T)  F).  .  =  lim  )  F  ( t  ,  x(l))  F  (t  ,  x(j)) 

-  -  Lj  L  LXU  U.  -  /  v  V*  -  ' 


(ID 


r\(t,  t')  F(t,  x(l))  F  ( t ' ,  x(j))  dt  dt' 


(FT  T)  S).  =  lim  Y  n  F ( t  ,  x(l))  S(t.  ) 

—  -X  -  l  Z_j  u,u  u  <■** 


(12) 


r 
I 

j  j 


n(t,  t')  f (t,  x(l))  s(t')  dt  dt7 


The  superscript  T  denotes  a  vector,  matrix,  or  operator  transpose. 
h(t,  t')  is  the  operator  inverse  of  cp ( t ,  t'),  here  written  symbol¬ 
ically  as  an  integral  kernel. 

Now  let  0  be  any  real-valued  functional  of  o/  having  a  finite 
second  moment  with  respect  to  the  a  priori  p.d.f.  of  o'-  Then  0,  the 
Bayes  optimum  estimate  of  0,  is  the  expected  value  of  0  with  respect 
to  p  (a  |  S) .  In  particular,  let  linear  and  quadratic  functionals  be 
defined : 


0.  =Y  a(x(l))  L(x (l) )  =  L  •  o'  (13) 

I  .  /  I  -M- 

i 

0n  =  Y  cy(x(i))  *(x(j))  Q(x(i),  x(j))  =  o,T  Q  a  (14) 

z_i  *»  —  ■**  —  — 

i  >  J 

Then, 
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§L  -  I  ^ 


})  a(x(l))  =  L  •  a 


(15) 


=  Trace 


[b-1  2 


+  a  Q  o? 


(16) 


These  formulas  are  obtained  simply  by  taking  the  expected  values 

r\ 

of  f3  or  0  with  respect  to  p(a  |  S)  .  Note  that  (3  is  not  formed 

a  C  i ) 

only  by  inserting  the  estimates  a(x  )  into  the  expression  for  0^; 
an  additional  term  is  necessary.  In  a  similar  manner  the  Bayes 
estimates  can  easily  be  written  down  for  third,  fourth,  or  higher 
order  forms  in  a,  applying  well-known  formulas  for  higher  moments  of 
Gaussian  distributions  in  terms  of  their  means  and  covariances.  These 
formulas  would  be  applied  to  find  expected  values  with  respect  to 
p(a  |  S);  the  resulting  estimates  would  be  functionals  of  a  and  B  ^ . 

These  constitute  a  class  of  nonlinear  estimation  cases  in  which  solutions 
valid  for  all  signal- to-noise  ratios  are  available. 


B.  DEPENDENCE  ON  ADDITIONAL  PARAMETERS 

It  is  of  interest  to  see  what  happens  when  the  signal  statistics 
depend  nonlinearly  on  an  additional  unknown  parameter  vector  X •  This 
may  happen  if  F(t,  x) ,  A_  ,  or  cp(t,  t^  depend  on  additional  para¬ 
meters.  We  will  be  chiefly  concerned  with  the  case  where  F(t,  x) 
depends  on  but  neither  nor  cp ( t ,  t;)  do  so.  The  approach  to  be 

taken  will  reflect  the  effect  of  such  additional  unknown  parameters 
on  the  Bayes  estimates  of  a  and  will  also  consider  the  question  of 
optimum  estimation  of  the  components  of  X,  itself. 
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In  general,  let 


P(S, 


a, 


X)  =  joint  a  priori  p.d.f.  of  S,  and  X 


(17) 


p(S,  a,  x)  =  p (S  |  a,  X)  p (o'  |  X)  q (X)  (18) 

Here,  S  is  considered  finite-dimensional.  The  function  q(X)  is  the 
a  priori  p.d.f.  of 
Now,  define 

p (a,  X  |  S)  =  joint  a  posteriori  p.d.f.  of  (a,  X)  given  £  (19) 

P(S,  a,  x) 

=  p— - 

!  I  P(S,  a,  X)  d&  dx 


p(a  |  S)  =  marginal  a  posteriori  p.d.f.  of  a  given  S 


(20) 


I  p(a,  X  S)  dX 

J  a**  1  '*** 


p(a  |  S,  X)  =  conditional  a  posteriori  p.d.f.  of  a  given  S  and  X  (21) 


P(S  |  a,  a)  p(a  |  X) 

[  P(S  |  a,  X)  p(o'  I  X)  da 


q(S  x)  =  conditional  p.d.f.  of  S  given  X 

1  AM  ’  MM 


f  p (S  I  a,  X)  p(a  I  X)  d£ 


(22) 
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q  (x 

|  S)  =  a  posteriori  marginal  p 

.  d .  f 

•  of  X 

given  S 

(23) 

t. 

■» 

p(a,  X 

1  §)  da 

q(s  ! 

x)  q(x) 

l 

1  q(S  | 

X)  q  (X)  d\ 

***  “• 

&(X) 

=  cond 

i  tional 

Bayes  estimate 

of  p 

given 

X 

(24) 

=  0(a)  p  (a  | 

S,  \)  da 

/W,  -n 

0  = 

overall 

Bayes  optimum  estimate 

of 

0  (cv) 

Am* 

(25) 

= 

r  ^  (o?)  i 
, .  ~*\ 

p(a,  X  1 

S )  da  dX 

XA, 

= 

[  P(o)  1 

o(a  1  S) 

"'"A  *** 

da 

Then,  it 

is  veri: 

fied  by 

direct  substitu 

tion 

that 

P  -  j  3(X)  q(X  I  S)  ciX  (26) 


j  &(X>  q(S  |  X)  q  (X )  dX 

» 

q(s  |  x)  q(x)  d\ 
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Equation  (26)  says  that  the  overall  Bayes  estimate  is  the 
expected  value  of  the  conditional  Bayes  estimate,  given  with 
respect  to  the  marginal  a  posteriori  p.d.f.  of  \  given  S.  However, 
the  function  of  q(S  I  X)  can  be  determined;  this  is  done,  in  fact, 
by  following  the  same  procedure  used  to  derive  Eq .  (5),  carefully 
keeping  track  of  all  the  terms  involved. 

The  result  is 


q(S  I  P  =  p(S  I  X)  p(a  |  X)  da 


_t 

2 


-N  ^ 

=  (det  A  B)  (2tt)  "2  (det  cp) 


x  exp 


{-* 


r  T  T  -IT' 

S  q  S  -  S  T}  F  B  F  q  S 


1 


where  A,  B,  cp,  q,  F  may  be  functions  of  X-  Equation  (27)  applies  to 
the  case  where  the  observation  times  are  restricted  to  a  set  of  N 
points  (tj- 

Let  us  now  assume  that  cp,  the  covariance  function  of  the  additive 
noise,  is  independent  of  \.  Then,  when  Eq .  (27)  is  substituted  into 

%  ~  "I? 

Eq.  (26),  (det  ^)  and  (2n)  2  drop  out,  so  it  is  possible  to  go  to  the 

limit  as  t  \  becomes  dense.  Also,  the  following  formula  may  be  used: 
t  u ) 


a(\)  =  B_1(X)  FT(x)  q  S 

_  —  **  V*  *•*  ^  w  -V 


(27) 


(28) 
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The  result  is 


P  = 


J  ^  [ 


det  A(X)  B(X) 


-h 


exp 


r%  ;T 


ft)  B  (X)  aft) 


q  (X )  dx 


det  A(x)  B  (X ) J  exP  [h  »T(X)  B(X)  g(X)"j  q  (X)  dx 


(29) 


In  some  special  cases  which  will  be  of  interest  in  later 
sections  it  is  found  that  not  only  A  but  also  B  is  independent  of 
In  such  cases 


C  A 


J  P(X)  exp  gT(x)  B  q(X)j  q  (X)  dx 


{3  = 


exp  \h  *T(X)  B  a(X)}  q(X)  dx 


(29a) 


The  expression  Eq.  (27)  for  q(S  I  \)  can  also  be  used  to  deter- 
mine  Bayes  optimum  estimates  for  functionals  of  Thus,  suppose 
y (X )  is  any  real-valued  functional  of  Then,  if  r)  is  independent 


of  X 


-h 


A  Y(x>  [det  A(x)  B(x)J  exp  h  aT(X)  B(X)  g(X) 


y  - 


q(x)  dx 


■h 


[det  A(x)  B(x)  exp  [%  aT(X)  B(x)  a(\) 


q  (x)  dx 


(30) 


If  A  and  B  are  independent  of  \ 


I  ^ 

exp 

h  a  (X)  B 

a(x) 
_ — _ 

q  (X  )  dX 

1 6XP 

h  a 1  (X)  B  ct(\ ) 

q(x)  dx 

M 

(30a) 
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A  ^ 

C.  ALTERNATIVE  METHODS  OF  EVALUATING  a  and  P 

A  /» 

An  alternative  method  for  evaluating  ty  and  p  for  various 
functionals  of  can  be  obtained  from  linear  optimum  signal  esti¬ 
mation  theory.  Starting  with  Eq .  (1),  suppose  the  goal  is  to  estimate 

<y(x^^),  given  S(t).  The  correlation  function  of  with  S(t) 

is 

E  fa(x(l))  S(t)]  =  )  A  F(t,  x(j))  (31) 

j 

The  autocovariance  function  of  S(t)  is 


* 

<P  (t. 


t') 


=  E 


S(t) 


(32) 


=  cp(t,  t')  +  y  Ajk  F(t,  x(j))  F (t 7 ,  X(k)) 

j.k 


Then,  the  Bayes  optimum  estimate  cyCx^^)  is^ 


O' 


(i(i))  =IAu  llm  l  V  F(V  -<J))  slt^ 


M-  3 


(33) 


where 


q*  =  p,,  element  of  matrix  inverse  to  cp  (t^,  t^) 


p  v 


(34) 


This  may  be  written 


*  T 

cp  =  cp  +  F  A  F 

•/A  ,  ..  ‘AN  V*  *»*' 


*  *  -1 

T)  =  (9  ) 


(35) 

(36) 


A  T  * 

a  =  A  F  r\  s 


(37) 
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In  this  form,  it  is  immediately  apparent  how  these  expressions 
may  be  generalized  when  the  sum  in  Eq .  (1)  is  replaced  by  an  integral. 
Specifically,  suppose 


and  let 


S(t) 


|  cy(x)  F  (t,  x)  dx  +  s  (t) 
X 


(la) 


A(x,  x7)  =  E  |jy(x)  a(x') 


(2a) 


Then,  Eqs .  (35) -(37)  still  hold,  but  with 


F  A  F  = 


F(t,  x)  A(x,  x7)  F(t7,  x7)  dx  dx 7 


(38) 


X  X 


T 

A  F  = 


A(x,  x7)  F(t,  x7)  dx7 


(39) 


X 

and  so  forth. 

Still  another  alternative  expression  for  a  can  be  stated  for 
the  cases  represented  by  Eq .  (1),  using  the  finite  sum  expression, 
when  A  is  a  diagonal  matrix.  Suppose  it  is  desired  to  estimate  a(x^). 

The  sum  over  j  '/  i  can  be  treated  as  additional  noise.  Thus,  the 
problem  as  given  by  Eq .  (1)  can  be  replaced  by  an  equivalent  one  in 
which 


S  ( t)  =  c*(x(l))  F  ( t ,  x(l))  +  7(t) 


(lb) 


where 


cpi(t,  t7)  =  E  |^e(t)  e(t7)J  =  cp(t,  t7)  +  ^  A_  F(t,  x^)  F(t7,  x'^) 


(40) 


jA 
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This  gives  rise  to  an  equivalent  problem  in  which  the  matrix  Ih 
is  a  scalar.  T'io  formula  for  a(x^)  becomes 


tv 


(x 


(i) 


)  = 


f  F ( t ,  S(l))  q,  (t,  t')  SCt')  dt  dt' 

u  JL 

(An)  1  +  J  j  F(t,  x^)  T)i(tl  t ')  F  (t ' ,  x(l))  dt  dt' 


(41) 


where  the  integrals  are  symbolic  representations  of  limits  of  finite 
sums  over  p,,u  of  the  kind  used  above. 

In  many  cases,  the  matrices  cp^  ate  virtually  independent  of  i 
(i.e.,  the  addition  of  terms  involving  i  to  the  sum  in  Eq .  (40) 
would  be  negligible);  in  such  cases,  the  result  would  be 

cp  «  cp  +  F  A  FT  (42) 


and 


T 

..  (i), 

a(x  )  -  ■ 


a(x(l))  +  (FT  q  F) 


li 


where  a  is  defined 


(O  .«)■ 


-1 


(43) 


(44) 


Both  the  expressions  for  a  given  here  and  those  given  in 

Section  II. A  are  useful.  Those  in  Section  II .A  are  most  useful  if 

the  scattering  body  can  be  represented  as  relatively  few  scatterers 

(mo  e  precisely,  if  there  are  relatively  few  positions  x^  which  are 

candidates  for  the  position  of  a  major  scattering  center);  or  if 
T  - 1 

(F  tj  F  +  A  )  is  for  other  reasons  relatively  easy  to  invert.  On 
the  other  hand,  the  expressions  in  the  present  subsection  are  useful 
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if  cp  E  A  F  is  relatively  easy  to  invert  (say,  if  the  total  time 
of  observation  is  very  short). 

Finally,  it  is  of  interest  to  discuss  alternative  evaluations  of 
and  the  manner  in  which  these  generalize  to  the  case  in  which  S(t) 
is  given  by  an  integral  expression  of  the  form  Eq .  (la). 

A 

First,  when  S(t)  is  given  by  a  finite  sum,  3^  is  given  by 
Eq .  (16).  Either  of  the  expressions  Eq .  (9)  or  (37)  can  be  used  for 

A 

a  (if  A  is  diagonal,  Eq .  (41)  can  also  be  used),  but  B  is  still 
given  by  Eq.  (10).  Also,  if  Q  depends  on  only  one  point  x^  and 
A  is  diagonal,  then  B  in  Eq .  (16)  can  be  replaced  by  a  scalar: 


~  1  T  T 

B.  =  - -  +  F.  ri  .  FT 

l  A.  .  l'li 

li 


(45) 


The  manner  in  which  the  expression  Eq .  (16)  for  generalizes 
to  the  cases  where  S(t)  is  given  by  the  integral  expression  Eq .  (la) 
is  not  completely  obvious.  Here,  is  defined 

-  [  a(x)  Q  (x,  a(x  )  dx  dx  ‘ 


(14a) 


A 

There  is  no  difficulty  in  generalizing  the  term  a  Q  a  in  Eq-  (16)-- 
this  is  simply  carried  over  with  a  given  by  Eq .  (37).  The  problem 
is  in  deciding  how  to  generalize  B  in  Eq •  (16). 

This  problem  can  be  approached  as  follows.  Comparing  two 
expressions  for  a ,  Eqs .  (9)  and  (37),  it  can  be  seen  that 


B 


-1 


T 

F  r| 


A  F 


* 

n 


(46) 
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Thus 


- 1  T  *  T  - 1 

B  =  A  F  T|  F  (F  T)  F) 


(47) 


In  this  form,  the  expression  generalizes  to  the  cases  represented  by 

*"  1 

the  integral  expression  Eq .  (la).  The  interpretation  of  B  in  such 
cases  is 


B  1  (x,  x;) 


A (x ,  y)  F ( t ,  y)  q  (t,  O 
X  F(t/,  z)  U(z,  x  )  dy  dz  dt  dt/ 


(48) 


where 


U  1  (x ,  x  ')  = 


■»  if 

F(t,  x)  T]  ( t ,  t^)  F  ( t  ,  x^)  dt  dt 


(49) 


These  are  all  symbolic  expressions  for  operators. 

A 

Expressions  for  3  when  3  is  a  higher  order  form  in  a  can 
similarly  be  developed  for  the  integral  case  represented  by  Eq .  (la). 

One  can  regard  B  ,  given  by  Eq .  (47),  as  the  covariance  function  of 

AM 

the  (Gaussian)  conditional  distribution  of  a  given  S;  and  a,  given 

A 

by  Eq  •  (37),  as  the  mean-value  function.  The  expressions  for  3 
for  higher  order  forms  are  then  obtained  using  the  well-known 
expressions  for  the  higher  order  moments  of  a  Gaussian  random  process 
in  terms  of  its  means  and  covariance  function. 


D.  FOMENTS  OF  ESTIMATION  ERRORS 

Often  the  problem  of  primary  interest  is  to  determine  the  mean 
square  estimation  error.  Of  course,  given  the  expressions  for  the 
estimates  themselves,  the  mean  square  error  can  be  estimated  by 
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Monte  Carlo  techniques.  It  is  much  preferable,  however,  to  use 
analytic  expressions. 

Such  expressions  can  be  stated  for  linear  and  quadratic  functionals 
of  a;  it  is  easy  to  see  how  they  could  also  be  stated  for  higher  order 
functionals  by  extending  the  same  techniques  in  an  obvious,  though 
tedious,  manner.  These  statements  are  made  assuming  r.o  dependence 
of  F  on  additional  unknown  parameters.  If  such  dependence  exists, 
the  problem  of  writing  expressions  for  mean  square  error  is  much 
more  complicated;  of  course,  lower  bounds  can  be  obtained  by  using 
the  expressions  which  assume  the  unknown  parameters  to  be  known 
(or  by  using  Cramer-Rao  inequalities). 

The  mean  square  error  expressions  of  interest  are  the  conditional 
mean  square  error,  given  that  a  has  a  definite  value  a  ;  and  the 
unconditional  mean  square  error,  given  an  arbitrary  "true"  statistical 
distribution  of  o’-  Cases  of  interest  must  include  those  in  which  the 
"true"  statistical  distribution  of  cv  differs  from  the  a  priori  distribu¬ 
tion  assumed  for  purposes  of  forming  the  Bayes  estimate,  although  a 
special  case  of  interest  arises  when  these  two  distributions  coincide. 

It  may  also  be  of  interest  to  determine  mean  square  error  when  the  true 
covariance  of  the  noise  [e(t)}  differs  from  that  assumed  in  forming 
the  Bayes  estimate. 


First  consider  the  means  and  second  moments 
conditional  on  a  =  a  •  recall  that 

a**  O 


of  PL  and  Pq 


S  =  F  a  +  e 


(50) 
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Thus,  using  the  formulas  of  Section  II  .A  for  a, 

E  I  PT  I  a  1  =  L  B_1  FT  t)  F  a  (31) 

L  L  -wO  AM  ^ 

and 

E  F  (P  -  PT)2  |  a  |  =  L  B_1  FT  t)  9  T]  F  B"1  LT  (52> 

L  L  t,  woj  ~  -  -  J  „o  -i  ~  ^ 

+  |  L  B_1  FT  ti  F  y  -  0  (Of  )1 

Here,  cpQ(t,  t7)  is  the  "actual"  covariance  of  e;  of  course,  if  this 
is  also  the  covariance  assumed  in  the  Bayes  estimation,  then 
T]  tpn  B  =  T]  • 

If  the  assumed  a  priori  distribution  of  a  has  very  large 

variances,  i.e.,  if  the  estimates  depend  mainly  on  the  observed 

signal  and  only  weakly  on  the  a  priori  mean  values,  then 
-1  T 

B  F  ti  F  a  I  and  the  second  term  in  Eq .  (52)  vanishes,  while  the 

*+  t+m,  M 

right  side  of  Eq  .  (51)  becomes  simply  j3  ( q *  ). 

A 

The  moments  of  P_  conditional  on  a  =  a  can  best  be  determined 
Q  —o 

as  follows.  First 

§n  -  Pn(a  )  =  eT  M  e  +  R  e  +  W  (53) 

Q  Q  -  —  ~  - 

where 


n  F  B"1 

-1  T 

Q  B  F  T) 

(54) 

««*  ■** 

T  T 
2  a  F 

M 

(55) 

*»o  *** 

•*v> 

T  „T  „ 
a  F  M 

Fa  -  Pn(a  ) 

(56) 

**  *■* 

**  ~o  Q  ~o 

r  -1 

+  Trace 

LB-  2J 
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Therefore, 


=  Trace 


-1  "1 

B  Q  +  a 

-  -  ~c 


+  Trace  i  M 


$ 

*o_ 


F  M  F  a 

****•>••*  v*Q 


and 


E  |^(eT  M  e)2  |  aQj 

+  E  L(E  i)2  !  2„] 


+  2  W  E 


M  e 


“1 


Using  standard  formulas  for  fourth  order  moments  of  Gaussian 
processes  , 


L(BQ 


"  PQ)2  1  ~ol  =  LTrace  i? 


+  2  Trace  (M  cp  M  cp  ) 
~  Jo  -  Zo 

-i-  R  cp  RT 
-  -o  - 

+  2  W  Trace  [m  cp  ] 


(57) 


(58) 


(59) 


Now,  suppose  it  is  desired  to  obtain  the  unconditional  expectations 

with  respect  to  some  given  "actual"  statistical  distribution  of  a  , 

~o 

which  can  be  denoted  dP(aQ)*  All  such  moments  are  simply  the 

integrals  of  the  conditional  expectations  with  respect  to  dP (o'  ). 

—  o 

If  dP(c)'  )  is  Gaussian  with  covariance  A  ,  then  the  procedure  to  be 

~o  -o 

followed  is  similar  to  that  for  determining  the  conditional  moments. 
Thus,  all  the  conditional  moments  above  can  be  expressed  as  the  sum  of 
constants  and  linear  and  quadratic  forms  in  a  •  The  moments  of  all 
these  can  be  expressed  in  terms  of  the  matrix  A  .  For  example 
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«Q  -  V2] 


Trace 


<2  2o>l 


+  2  Trace  (M  cp  M  cp  ) 
-  -o  -  -o 


-1 


+  2  Trace (M  9  )  Trace (B  Q 
—  -o  ~  •* 

+  4  Trace  f*A  FT  M  cp  M  F 
L~o  -  -  -o - J 

r  T 

+  Trace  j  A  (F  M  F  -  Q) 

L“  O'"-  ^  ^  _ 


(60) 


These  expressions  all  can  be  generalized  to  the  case  where  a  is  given  by 
Eq .  (37),  and  hence  to  the  case  where  S(t)  is  given  by  the  integral 
Eq .  (la),  simply  by  interpreting  B  ^  in  all  the  above  expressions 
via  Eq .  (47). 


These  techniques  can  obviously  be  extended  to  determine 
the  characteristic  functions  (and  hence  the  p.d.f.'s)  of  the 

A  . 


estimates  3T  and  p^,  as  well  as 
order  forms  in  a.  However,  the 
tedious  . 


the  moments  of  estimates  of  higher 
necessary  algebra  becomes  extremely 


E .  ADDITIONAL  REMARKS 


It  is  clear  how  the  function  Q(x,  x  )  can  be  chosen  so  as  to 
yield  information  about  the  integrated  radar  cross  section  contributed 
by  specified  sections  of  the  target,  as  well  as  to  make  estimates  of 
body  dimensions  (e.g.,  by  letting  Q(x,  x  )  31  x  6  (x  -  x),  where 
a  is  a  large  integer).  It  will  also  be  of  interest  in  this  connec¬ 
tion  to  form  estimates  of  other  functions  of  P  ,  or  of  functions 
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depending  on  two  quadratic  functionals.  For  example,  functions  of 

,1/m 


the  form  I  3 


Lpq] 

of  interest. 


or  3n  /3n  or,  more 

W1 


generally,  f  3  ,  3  may 

L  h!1  W2J 


be 


In  many  such  cases,  even  though  a  strong  signal  assumption  may 
not  be  justified  m  determining  from  o',  such  an  assumption 

.  /  ^  A 

may  be  justified  for  determining  f  by  f  =  f (  3„  ,  3n  )■  An  error 

41,  ^2 

/  A  A 

analysis  can  then  be  carried  out  by  expanding  f(  3  ,  3  j  to  first 

^1  ^2 

A  A 

order  in  the  error  terms  3„  -  3n  ,  3n  -  3  and  then  applying  the  error 

^1  Q1  g2  Q2 

A  A  A 

formulas  for  3n  and  3„  •  (This  would  require  cross-moments  of  3 

^1  Q2  kl 

A 

and  3„  ,  which  can  be  determined  by  the  same  sort  of  technique  as 


used  in  Section  II. D.) 

A  possibly  useful  method  of  evaluating  the  expressions  for 

A  i*  A 

a,  3L,  3q,  and  so  forth,  alternative  to  those  stated  above,  is  the 
use  of  recursive  Bayes  estimates . ^  The  application  of  recursive 
methods  to  the  problems  considered  in  this  Memorandum  is  currently 
under  investigation. 
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III.  APPLICATION  TO  ROTATING  TARGETS 

A.  TARGET  MODELS  AND  SIGNAL  MODELS 

The  target  and  signal  models  to  be  considered  can  be  represented 
in  the  following  way  (omitting  the  noise  term  e(t)) 

V  r  "1 

S(t)=;a(r.,cp.)cos|u)  t  +  r  sin(tu  t  +  cp. )  +  £  . 

Lx  i  L  c  i  r  Ti  lj 

i 

where 

UJ  =  carrier  frequency  (radians/ sec ) 

u)  =  target  rotation  frequency  (radians/sec) 

5 .  =  RF  phase  parameter  of  signal  from  itn  scatterer 
(further  discussed  below) 

(r.  ,  cp.)  =  location  of  i*”^  scatterer  in  a  coordinate  system 
i  i 

fixed  to  the  rotating  body  (further  described 
below) 

a(ri,  cp  )  =  amplitude  (absolute  magnitude)  of  i ^  scatterer 
In  what  follows,  cp  always  represents  an  angular  coordinate  unless 
otherwise  specifically  stated  and  should  not  be  confused  with  the 
covariance  function  cp  of  the  observation  noise. 

The  coordinate  system  (r,  cp)  is  defined  as  follows.  Initially, 
suppose  that  the  axis  of  rotation  is  normal  to  the  line  of  sight 
to  the  radar.  Then,  the  (r,  cp) -plane  is  normal  to  the  rotation 
axis  and  contains  the  line  of  sight.  The  coordinates  (r,  cp)  are  then 
cylindrical  coordinates  of  a  point  on  the  target  in  a  coordinate 
system  fixed  to  the  body,  in  which  the  axis  of  rotation  is  the  z-axis. 
At  any  time  t,  the  angle  of  the  9=0  axis  with  respect  to  a  plane 
normal  to  the  line  of  sight  is  (A)  t.  (In  other  words,  oj^  t  +  n/2  is 
the  angle  between  the  cp  =  0  axis  and  the  line  of  sight.)  The  units 


(61) 
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of  r  are  chosen  so  that  one  wavelength  equals  4tt  units.  Otherwise 

a  factor  4tt/\  would  have  to  multiply  r  in  Eq .  (61). 

The  same  model  can  be  applied  to  some  types  of  targets  even  when 
the  axis  of  rotation  is  not  normal  to  the  line  of  sight.  Specifically, 

suppose  all  the  significant  scatterers  on  the  body  are  known  a  priori 

to  lie  in  a  plane  which  is  normal  to  the  rotation  axis.  In  such  a 
case,  let  (r,  cp)  still  be  cylindrical  coordinates  fixed  to  the  body, 
with  the  z-axis  being  the  rotation  axis.  The  signal  can  still  be 
described  by  Eq .  (61)  if  r^  is  replaced  by  r^  sin  Y,  where  Y  is  the 
angle  between  the  rotation  axis  and  the  line  of  sight.  In 
this  case,  t  is  the  angle  at  time  t  between  the  cp  =  0  axis  and 

the  line  of  intersection  of  the  plane  normal  to  the  line  of  sight 
and  the  (r,  cp)  plane  fixed  to  the  target.) 

The  symbol  ^  denotes  an  additional  unknown  phase  parameter 

th 

associated  with  the  signal  reflected  from  the  i  scatterer.  It  is 
assumed  to  be  uniformly  distributed  over  (0,  2n) •  Thus,  Eq .  (61)  can 
also  be  written 


S  ( t)  =£a(r.,  cp^ ,  1)  sin  t  +  r_^  sin(cr^  t  +  qO 

i 

+  y  a(r.,  q> .  ,  2)  cos  I  uj  t  +  r.  sin  (u)  t  +  cp, ) 

»  1  1  L  c  l  r  1  _ 


where 

C«'(ri,  q^,  1)  =  a(ri,  qh)  cos  ^ 
O' (r± ,  cp.,  2)  =  a(ri?  cpi)  sin  ? 


(62) 


(63) 
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Thus,  the  superscript  i  of  Section  II  is  here  a  pair  of  indices 

(i,  (i)  where  i  indexes  the  scatterers  and  p  takes  on  only  two  values, 

1  and  2.  Also 


x  =  (r,  cp,  p)  (64) 

where  p  =  1  or  2 .  It  is  assumed  that  the  components  of  q.  have  a  priori  a 
joint  Gaussian  distribution. 

Since  there  are  no  modulating  functions  multiplying  the  sum  in 

Eq .  (61),  CW  illumination  of  the  target  is  in  effect  being  assumed. 

However,  the  same  model  can  also  represent  the.  case  of  a  pulse  radar 
for  which  the  range  resolution  of  a  pulse  is  much  larger  than  the 
physical  extent  of  the  target,  and  for  which  the  pulse  repetition 
frequency  is  sufficiently  high  that  the  maximum  distance  moved 
by  any  scatterer  in  the  interpulse  period  is  small  compared  with  a 
wavelength . 

It  will  be  assumed  that  yj^  is  known  a  priori.  This  is  not  an 
unrealistic  assumption  in  many  cases  (e.g.,  when  the  target  is  known 
a  priori  to  be  stabilized  and  the  effective  rotation  is  due  to  non- 
radial  translatory  motion  of  the  target  with  respect  to  the  radar); 
also,  it  can  often  be  assumed  that  the  rotation  rate  has  bean 
accurately  established  from  prior  observations,  since  rotation  rate 
is  a  parameter  often  capable  of  being  very  accurately  established. 

It  is  clear  that  the  above  model  implies  that  a  coherent  radar 
is  used.  This  does  not  imply  that  signal  phase  is  known  a  priori,  but 
merely  that  the  RF  phase  of  the  received  signal  plus  noise  is  observed 
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and  made  use  of.  An  interesting  problem  is  to  determine  the  effect 

of  lack  of  coherence,  either  in  the  radar  processing  or  because  of 

phase  errors.  It  is  known  that  much  signature  analysis  can  be  done 

even  with  completely  incoherent  operation  (e.g.,  even  if  the  radar 

measures  only  signal  amplitude).  The  discussion  of  this  question 

is  given  in  Section  III .E  below. 

The  targets  envisaged  by  the  representation  stated  in  Eq.  (61) 

consist  of  sets  of  rigidly  connected  point  scatterers.  The  potential 

positions  of  the  point  scatterers  with  respect  to  the  rotation  axis 

are  assumed  to  be  known  a  priori  to  within  approximately  one-half 

wavelength  (one  wavelength  in  the  two-way  path).  The  additional  unknown 

til 

parameters  ^  are  included  to  allow  the  precise  position  of  the  i 

scatterer  to  differ  from  the  position  (r^ ,  q_^)  by  up  to  one-half 

wavelength.  (For  some  further  cases  treated  below,  it  is  required  that 

all  the  ^'s  be  equal  or  highly  correlated.)  It  is  necessary  to  assume 

that  the  maximum  separation  between  potential  positions  (r^,  cpu)  of 

scatterers  is  of  the  order  of  one-half  wavelength.  Strictly  speaking, 

th 

(r\  ,  qh)  should  be  called  the  i  mean  potential  scatterer  position. 

It  is  somewhat  easier  to  visualize  the  potential  scatterer 
positions,  projected  on  the  (r,  cp)-plane,  as  a  regular  rectangular 
lattice,  with 

.  th  .  .  ,  s 

1  position  =  (u  ,  v  )  (63) 

r  m  n 

i  -  (m,  n) 

u  =  r.  cos  cp. 

m  l  l 

v  =  r.  sin  cp. 

n  i  i 
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On  occasion  the  corresponding  integral  form  over  the  (r,  cp)-plane 
will  be  used. 

The  scattering  model  just  stated  also  applies  to  sets  of  rigidly 
connected  spheres,  even  when  the  spherical  radii  are  larger  than  a 

Vc 

wavelength.  In  such  cases,  the  scattering  point  on  each  sphere  is, 

at  any  instant,  the  point  on  the  sphere's  surface  intersected  by  the 

line  connecting  the  center  of  the  sphere  to  the  radar.  If  the  center 

of  a  given  sphere  describes  a  circle  of  radius  r,  so  also  does  the 

instantaneous  scattering  point,  but  the  circle  is  displaced  toward 

the  radar  by  the  radius  of  the  sphere.  If  the  scattering  spheres 

have  different  radii,  this  is  accounted  for  by  the  phase  variables 

In  the  models  to  be  treated,  it  is  assumed  that  each  scatterer 

is  visible  during  the  whole  observation  interval.  This  might  be 

true  if  the  target  scatterers  were  all  located  on  a  rotating  disk,  or 

on  a  rotating  line,  with  slightly  tilted  rotation  axis;  or  if  the 

target  otherwise  consisted  of  individual  points  which  did  not  shadow 

each  other.  A  similar  model  might  apply  to  long  thin  bodies  which 

are  sufficiently  symmetrical  about  the  lengthwise  axis. 

In  many  cases,  one  is  interested  in  bodies  for  which,  even  if 

they  can  be  reasonably  modeled  as  collections  of  rigidly  connected 

point  scatterers,  each  scatterer  is  visible  only  part  of  the  time,  i. 

is  visible  over  some  restricted  range  of  values  of  cp  +  oj  t,  as  on 

r 

the  surface  of  a  rotating  sphere.  It  will  be  apparent  that  the  treat 
ment  to  be  given  can  be  generalized  to  such  cases,  but  detailed 
evaluations  are  much  harder  and  will  only  be  indicated. 

*It  should  be  remembered  that  (r^,  cfu)  are  possible  scatterer 

locations;  in  any  actual  case,  many  of  these  possible  locations  might 
be  empty,  that  is,  occupied  by  zero  amplitude  scatterers. 
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In  all  cases,  it  will  be  assumed  that  the  observation  noise  is 


cp(t,  t7)  =  6(t  -  t7) 

n (t,  t7)  =  ~  6 (t  -  t ') 
o 

Two  categories  of  target  models  will  be  treated:  first,  a  model 
in  which  the  potential  scatterer  positions  are  not  assumed  a  priori 
to  be  collinear;  and  second,  a  model  in  which  it  is  assumed  that  the 
scatterer  positions  are  known  a  priori  to  be  collinear. 


(66) 

(67) 


B.  NONCOLLINEAR  A  PRIORI  SCATTERER  DISTRIBUTIONS 

For  noncollinear  a  priori  scatterer  distributions,  the  matrix  A 
will  be  taken  to  be  diagonal: 


A(ri, 


Mm 


M> 


')  =  A, 


6ju 

ip 


(68) 


where  6  is  the  Kronecker  delta.  The  expressions  necessary  to  evaluate 
the  Bayes  estimates  can  then  be  derived  by  either  of  two  approaches. 
First  Approach 

The  first  approach  will  be  to  calculate  B(x,  x7)  by  application 
of  the  formula  Eq .  (7).  This  requires  evaluation  of  the  following 
types  of  integrals: 


11 


12 


sin  t  +  r  sin  (u^  t  +  ,)]  sin  [„c  t  +  r  sin  (u>^  t  +  cp  )]  dt  (69) 
sin  |^uc  t  +  r  sin  (u)r  t  +  cp)]  cos  u>c  t  +  r7  sin  (a)  t  +  cp")  dt  (70) 


‘22  J 


cos 


UJ  t  +  r  sin  (u)  t  +  cp) 
.  c  r 


[■ 


cos  |  t  +  r  sin  (0)^  t  +  cp7)”!  dt  (71) 
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=  same  as  with  r, 


cp  and  r  ,  cp  interchanged 


(72) 


The  integrals  can  be  assumed  to  extend  over  an  interval 
J 0  |  =  t|  =  8^;  without  loss  of  t°nerality  it  can  be  assumed  that 

t  =  0  is  the  midpoint  of  the  observation  interval. 

If  is  sufficiently  large,  as  it  will  be  in  virtually  all 
cases  of  interest,  then 


himl 


22  2u)  .! 

r 


cos  |^|  p  |  sin  (9  +  9  ) 


d6 


_  1  r 


12  2oj 


21  2co 


sin 


sm 


*  - 


p|  sin  (8  +  9  )  dB 


p  sin  (8  -  9  )  |  d8 


where 


P  = 


2/2  /  .  / . 
r  +  r  -  2r  r  cos  (cp  -  cp  ) 


*  -1 
0  =  tan 


.  /  .  / 
r  sm  cp  -  r  sm  cp 

<  r  cos  cp  -  r '  cos  cp7 


Now,  the  integrands  can  be  expanded 


cos 


p  sin  (0+8 


-  J  P  +  2 


>] 

00 

I  J2n  (lp!  cos 


n=l 


s  in 


*  n 


p|  sin  (0  +  0  ) 


-  2  )  J 


2n+l 


n=0 


r  *  - 

2n(8  +  0  ) 
p|^sin  j^(2n  +  1)  (8  +  8  )J 


Consequently,  recalling  that  the  integrals  in  9  =  id  t  are 

assumed  to  extend  from  8  =  -  9W  to  8  =  +  9,,, 

M  M 


(73) 

(74) 

(75) 


(76) 

(77) 


(78) 


(79) 


34 


11 


_  /  i 

*22  uj  J0  l|P 

r 


(JO 


I  h„0>0 


n=l 


cos  2n  9  sin  2n  9 


M 


n 


'12  -I21 


=-  y 

uj  Li 


J2n+1  HO 


n=0 


* 


sin  (2n  +  1)  9  sin  (2n  +  1)  9 


M 


2n  +  1 


(80) 


(81) 


If  9^  =  K-n,  where  K  is  an  integer,  which  is  to  say,  if  the 
observation  interval  is  an  integral  number  of  periods  of  the  rotation, 
then 


hi  =I22  Jo 

r 


*12  =  '^l  =  ° 


(82) 

(83) 


Ktt 

If  9^  =  — ,  with  K  an  integer,  then 


I11  I22  2m  Jo  ^ 1 P 


(84) 


X  =  -I  =  — 
x  12  21  oj 


I 

n=0 


1=12 


n 


2n  +  1 


s  in 


(2n  +  1)  6*]  J2n+1  (|plj 


(85) 


The  formula  for  B(x,  x  )  in  terms  of  1^,  etc.,  is 


B  (r ,  cp,  1,  r ' ,  cp',  1)  =  B(r,  cp,  2,  r',  cp / ,  2)  (86) 

r  -,"1  2  /  / 

=  |^A(r,  cp) J  +  —  I11(r,  cp,  r,  cp) ,  if  (r,  cp  =  (r  ,  cp  ) 

o 

=  I11(r>  ^ >  v' >  cp/)>  if  (r»  CP)  f  (r7.  cp7) 

o 


B(r,  cp,  1,  r ' ,  cp',  2)  =  I12(r,  cp,  r\  cp') 

o 


(87) 
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B(r,  9,  2,  r  ,  cp',  1)  =|~I21(r,  cp,  r ' ,  cp')  (88) 

o 


It  is  understood  that  r  and  cp  take  on  only  a  discrete  set  of 


values  r . ,  cp.  . 

l  l 


Second  Approach 

The  second  method  will  be  applicable  to  evaluation  of  Eqs.  (35)~ 

T 

(37).  The  basic  step  is  the  evaluation  of  F  A  F.  It  will  also  be 

assumed  that  the  density  of  possible  scatterer  locations  is  sufficient 

T 

to  enable  the  sums  in  the  expression  for  F  A  F  to  be  replaced  approxi- 
mately  by  integrals.  Thus 


T 

F  A  F 


{  P(r,  cp)  | 


cos  t  +  r  sin  (u.'^  t  +  cp) 


cos 


0) 


t '  +  r  sin  (u)^  t  +  cp)  J 


(89) 


!  t  ft 

+  sin  t  +  r  sin  (uj^  t  +  cp)  I  sin  ^  t  +  r  sin  (m^  t  +  cp)  ^  r  dr  dcp 


where 


tit  >ii 

A (r ,  cp,  p,,  r  ,  cp  ,  p,  )  =  0,  (r,  cp,  p,)  ■/  (r  cp  P  ) 


(90) 


=  P(r,  cp)  r  d  r  d  cp, 


cc.  u)  =  (r', 


It  is  understood  that  r  dr  dcp  may  actually  represent  a  finite  increment 
of  area,  equal  in  size  to  the  area  over  which  the  a  priori 
distribution  of  a(r,  cp)  is  significantly  correlated.  The  minimum 
size  of  r  dr  dcp  would  be  about  a  half -wavelength  square. 

Equation  (89)  may  be  rewritten 

/ 

u)(t-t)+rsin(cu  t+cp) 
c  r  T 

-  r  sin  (0)  t /  +  cp) 
r 


T 

F  A  F 


P(r,  cp)  cos 


(91) 


r  dr  dcp 
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If  it  is  assumed  that 


P(r,  cp)  =  P  (r ) ,  0  g  cp  <  2tt 


(92) 


then 


T 

F  A  F 


2n  cos 


uj  (t  -  t  )  P  (r )  J 


L  C 


J, 


rUtl  -  cos  (t  -  t ;)]  [ r  dr 


(93) 


P(r,  cp)  reflects  the  a  priori  knowledge  of  the  distribution  of 

T 

scatterer  intensities.  Some  explicit  expressions  for  F  A  F  are 

2 


n 


(1)  P (r )  =  P  r  exp 

o 


L-fc)J 


T  .  ^  n+2 

F  A  F  =  nr  P  cos 

**V>  0  0 


X  ,F 


n  +  2 


11  2 


[“c(t  -  c'>]  r  (f  +  1 

2 

r  r 

.  i-  .  _2_  [ 

’  ’  2  L 


1  -  COS  IU  (t  -  t  ) 

r 


(94) 


where  ^F^  is  the  confluent  hypergeometric  function.  For  n  =  0, 


T 

F  A  F 


tt  r  P  cos 
o  o 


<Dc(t  -  t  )J  exp 


-  r 


exp 


tL 


COS  cu  (t  -  t  ) 
r 


(95) 


(2)  P(r)  =  Pq  exp 


T 

F  A  F 


2tt  r  P  cos 
o  o 


3/2 


t') 


1  +  2r 


1  -  CO  SO)  (t  -  t;) 
r 


(96) 
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(3)  P(r)=P,0  =  r  =  r  and  P(r)  =  0,  r  >  r 

o  o  o 


F  A  F  =  2n  r  P  cos 
~  —  ~~  o  o 


x  -^2  j^l  -  cos  u)^(t  -  t;)  | 


[^c(t  -  t')J  J1  jro  \j  2  [1  -  cos  uJr(t 

-% 


-  t 


ttA 


)];  (97) 


The  result  just  stated  for  P(r)  =  Pq  exp 
further  interpreted  as  follows. 


Xo) 


can  be 


Deno  ting  t  -  t '  by  T, 


(98) 


(99) 


Thur 


exp 


~T~  COS  U)  T 

2  r 


* 

6  (T) 


L  r  N  f 

where  6  (t)  represents  a  periodic  Dirac  delta  function  with  period 

2tt.  The  approximation  indicated  in  Eq .  (100)  is  valid  in  the 

*  , 

following  sense.  If  both  the  6  function  and  the  function  on  the 
left  side  of  Eq .  (100)  are  regarded  as  kernels  of  integral  operators, 
the  result  ot  applying  these  two  operators  will  be  approximately  the 
same,  provided  the  function  to  which  they  are  applied  has  a  Fourier 

expansion  (over  toe  interval  2-rrA^)  which  is  negligible  for  all 

.  v  2 

n  *  n  where  n  <  <  H  r  • 
o  c  o 


(100) 
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Estimates  a(r^>  cp^)  involve  application  of  this  operator  via 
Eq .  (37)  to  F(t,  r  ,  q;  )  .  This  function  in  turn  satisfies  the 
criterion  just  stated  only  if  r.  <  <  r  .  Thus,  use  of  expression 
Eq .  (100)  would  be  justified  for  estimates  of  a  and  functionals  of 

v* 

a  only  for  a  region  r  <  <  r  . 


C.  COLLINEAR  A  PRIORI  SCATTERER  DISTRIBUTIONS  WITH  NO  SPECULAR  FLASH 
The  target  model  in  the  cases  now  to  be  considered  is  given  by 


S(t)  =  /  a(r.)  cos  1'^  t  +  r .  sin  (a)  t  +  8  )  +7,1 
L  l  L  c  l  r  o  l 

i 

or  equivalently 


(101) 


S(t)  =  )  cKr.,  1)  sin  |  a)  t  +  r.  sin  (uj  t  +  0) 

U  l  t-c  l  r  oj 

i 

r1  r 

+  ^  a(r^,  2)  cos  t  +  r^  sin  (u>^  t  +  8^)  | 
i 


(102) 


This  is  the  same  as  in  Eqs .  (61)  or  (62),  but  with  cp.  =  0.  Note  that 
r^  can  take  on  both  positive  and  negative  values  (otherwise  it  would 
be  necessary  to  admit  two  values  of  cp^)  • 

The  phase  factors  are  still  assumed  to  be  random  and  independent 

for  different  indices  i;  this  is  what  is  meant  by  the  statement  that 
no  specular  flash  is  assumed,  a  priori,  to  occur.  This  is,  however, 
an  assumption  about  the  state  of  a  priori  knowledge,  and  does  not 
preclude  cases  in  which  a  specular  flash  may  occur  in  actual  fact. 

A  specular  flash  is,  in  this  terminology,  said  to  occur  if  there  is 
some  viewing  aspect  at  which  all  the  scatterers,  or  a  dominant 
portion  of  them,  are  in  phase. 
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The  parameter  9^  is  the  target  aspect  at  t  =  0  and  is  considered 
to  be  an  a  priori  unknown  parameter.  In  the  notation  of  Section  II. B, 
•v  =  0  •  A  is  still  assumed  to  be  given  by  Eq .  (68). 

The  first  approach  of  Section  III.B  can  be  followed  with  the 
result 


B(r,  1,  r',  1  |  6o)  =  B (r ,  2,  r',  2  |  0q) 


=  [A (r ) ] ” 1  +  In(r,  r  |  0o),  if  r  =  r' 
o 


(103) 


■S'  hi(r-  r'  I  V  ifr/r 

o 


B  (r ,  1,  r  ,  2  |  0q)  =  -B(r,  2,  r' ,  1  j  0q) 


(104) 


r1!)*1'  I  V 

o 


where  1^  and  I  ^  ate  given  by  Eqs .  (80)  and  (81)  with 


/  i 

r  -  r  j 


* 

e  =  e 


(105) 

(106) 


In  particular,  if  0^  =  Kn,  i.e.,  if  the  observation  interval  is  an 
integral  number  of  revolutions,  then 


ll(r'  r'  ! 

e  )  =  —  J  (I  r  -  r'  1) 

o  UJ  o'  1 

r 

(107) 

12(r-  r' 

o 

II 

o 

CD 

(108) 
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Thus,  in  the  case  0..  =  Ktt,  B  does  not  depend  on  9  .  Therefore 

M  o 

Eqs .  (29a)  and  (30a)  of  Section  II. B  are  applicable  in  these  cases. 

9^  has  the  same  meaning  as  in  Section  III.B. 

Finally,  when  0,,  =  Kn 
M 


O'T(0  )  B  a(9  ) 
—  o  -  -  o 


S(t)  S(t/)  cos  I  id  (t  -  l  )  +  r.  sin  (i  t  +  9  ) 
L  c  i  r  o 


-  r.  sin  (uj  t;  +  9  )  dt  dt^ 
J  r  o 


This  defines  the  a  posteriori  p.d.f.  of  9  given  S  via  Eq .  (29). 

D.  COLLI NEAR  A  PRIORI  SCATTERER  DISTRIBUTION  WITH  A  SPECULAR  FLASH 

A  specular  flash  would  occur  at  some  target  aspect  if,  for  example 

^  5  for  all  i.  However,  xf  this  assumption  were  made,  the  joint 

probability  distribution  of  ^(r^,  1)  and  Q'tr^,  2)  would  not  be 

Gaussian  (even  though  the  marginal  distributions  for  any  one  index 

i  would  be)  unless  also  a^  =  a  P(i),  all  i,  where  P(i)  is  an  a  priori 

known  function.  Here  a,  has  been  written  for  a(r.). 

l  i 

The  latter  is  not  an  altogether  unrealistic  requirement,  since 
one  would  generally  expect  that  the  existence  of  an  aspect  where  the 
scatterers  are  all  in  phase  would  not  be  accompanied  by  a  very 
irregular  variation  of  amplitudes  of  the  individual  scatterers. 

For  example,  the  field  from  a  uniformly  illuminated  antenna  would 
be  equivalent  to  a^  =  constant. 

This  suggests  a  model  in  which  the  signal  from  a  linear  object 
giving  a  specular  flash  is  represented  by  Eq .  (101),  but  with 


(109) 
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l±  =  5,  all  i  (110) 

where  Z  now  represents  the  phase  of  a  signal  reflected  from  r  =  0, 
and 

ai  =  a  P (i ,  \  )  (111) 

* 

Wl‘erG  ^  is  a  vector  of  unknown  parameters  having  a  small  number  of 
components,  such  as  two  or  three. 

In  such  a  case,  o  is  a  two-component  vector  with  a  Gaussian 
a  priori  p.d.f. 


=  a  cos  5 


o2  =  a  sin  EJ 


A  can  be  assumed  to  be  a  diagonal  matrix,  with 


A11  -  A22  “  *  E<a  > 


A12  A21  ° 

i<  i< 

The  2x2  matrix  B(0o>  \  ),  with  components  B„(0o>  X  )  ,  is  derived 

via  the  first  approach  in  Section  III.Bo  To  illustrate,  consider 

the  case  8W  =  Krt,  where  the  observation  is  an  integral  number  of 
M 


(112) 


(113) 


periods.  The  result  is  then 


B  (0  ,  X*)  =  B(X*) 

**  o  *-> 


(114) 


?Q> 


2  ,  Ktt 


E  (a 


27  +  -  l  P<V  i  >  P(rj’  t  >  Vl  I) 

)  r  .  .  J  J 


(115) 


B12(X*)  =  B21(>.*)  =  0 


(116) 


Also,  if  \  =  (0o,  X  ), 


Q'1  (X)  =  Q'1  (9o,  X  ) 


(117) 


Bn(-*) 


-1 


r*  *  r  n 

,  P(r.,  X  )  sin  I  UJ  t  +  r.  sin  (uj  t  +  6  )  I  S(t) 
Z_i  i  ***  Lc  i  r  o 


dt 


(}s)  =  ^2(0o’ 


(118) 


,2  \ 


B. 


N  I  L  22 

O' 


a*)l  I  )  P(r.,  X  )  cos  j  to  t  +  r.  sin  (uu  t  +  0  )  S(t)  dt 
—  l-  Lc  i  r  o 


(P  5L>  2<P  LBllX*)]  i  j  S(t)  S(t')  < 1  )  P(r  X*)  P(r.,  \*) 

]\J  -I  ,J  ‘J  J  L—i  1  ***  J  AM 


J 


(119) 


cos 


UJ  (t  -  t')  +  r  sin  (iu  t  +  0  )  -  r.  sin  (uu  t '  +  0  ) 
L<-  1  r  oivr  o 


dt  dt ' 


It  would  be  necessary  to  use  Eq .  (29)  or  Eq.  (30)  in  this  case,  but 

the  expression  for  det  j  B (X )  is  very  simple: 

Am  J 


det 


B  (X) 


2 

r  *  -i 

=  Lbii(^  >J 


(120) 
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Since  the  total  number  of  unknown  parameters  in  this  case  has 
been  reduced,  by  assumption,  to  a  rather  small  number--^,  ,  6q,  and 
the  components  of  \--the  USe  of  the  infoi  ation  matrix  to  compute 
lower  bounds  on  estimation  errors,  and  the  use  of  maximum  likelihood 
estimates,  may  be  an  adequate  treatment. 

It  would  also  be  possible  to  consider  mixed  cases,  i.e., 
scattering  bodies  which,  a  priori,  have  a  component  of  the  type 
considered  in  this  subsection,  with  a  specular  flash;  and  another 
component  of  the  type  treated  in  Section  III.B.  It  would  be 
possible  to  apply  a  mixture  of  the  approaches  described  above.  For 
example,  for  purposes  of  estimating  \  or  0^  for  the  component  which 
provides  the  specular  flash,  it  would  be  possible  in  some  cases  to 
consider  the  other  component  to  provide  the  equivalent  of  noise 

■k 

(changing  the  noise  inverse  covariance  matrix  from  r)  to  ij  ,  but 
integrating  in  Eq .  (35)  only  over  the  scatterers  in  the  second 
component).  In  certain  cases,  the  resulting  equivalent  noise  could 
still  be  considered  white;  see  Eqs .  (98)-(100). 

E.  PHASE  ERRORS  AND  INCOHERENT  PROCESSING 

There  are  two  possible  sources  of  reduction  in  the  extent  of 
coherence  maintained  over  the  observation  time:  first,  phase  errors 
(such  as  those  due  to  propagation  medium  inhomogeneities,  or  to  phase 
jitter  in  the  radar,  or  to  uncorrected  platform  motion);  and  second, 
failure  of  the  radar  to  attempt  to  preserve  phase  coherence  over  the 
observation  time.  An  extreme  case  of  the  second  type  arises  when 
the  radar  observes  only  the  amplitude  of  the  returned  signal. 
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It  is  well  known  that  considerable  signature  analysis  is  possible 
even  with  amplitude-only  observations,  and  it  is  of  interest  to  com¬ 
pare  the  results  achieved  with  those  achieved  by  fully  coherent 
processing.  The  effect  of  phase  errors  on  coherent  processing  is 
also  of  interest.  The  following  discussion  describes  a  method  of 
approach  to  the  analysis  of  these  questions;  however,  many  questions 
remain  open,  even  with  regard  to  the  methodology. 

Let  us,  for  the  sake  of  argument,  revert  tc  the  signal  model  of 
Sections  III. A  and  III.B,  i.e.,  to  the  a  priori  noncollinear  collection 
of  scatterers.  However,  the  signal  (free  of  additive  noise)  will  now 
be  represented  as 


S(t)  =  a  (r^,  cpj  cos  jjJJ^  t  +  r^  sin  (u)^_  t  +  qO  -L  +  Q  (t)  J 


=  £  a(r.  ,  cp.  ,  1)  sin  jjj^  t  +  r±  sin  (u^  t  +  cp±)  +  C(t) 


(121) 


+  )  o(r.,  cPj  ,  2)  cos  |  iu  t  +  r.  sin  (uo  t  +  cp,  )  +  C (t)  I 
Z-j  ii  Lc  i  r  l  J 

where  (([(t)}  is  an  a  priori  unknown  random  process. 

The  distinction  between  the  phase  factors  and  (£(t)}  is  as 
follows.  The  quantities  E  are  fixed  for  all  time  and  represent 
relative  phases  of  the  (rigidly  connected)  individual  scatterers. 
The  quantities  £(t)  are  the  same  for  all  i,  and  represent  phase 
deviations  common  to  the  whole  scattering  body,  such  as  those  that 
would  be  produced  by  motions  along  the  line  of  sight  of  the  center 
of  rotation  of  the  body.  This  formulation  does  not  eliminate 
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the  contribution  of  the  phase  terms  r.  sin  (jj  t  +  cp.  )  to  the  informa- 
tion  which  can  be  extracted  from  signature  analysis,  even  if  the 
process  LG(t)}  is  uniformly  distributed  and  has  very  short  correlation 
time.  (If  the  phases  were  represented  as  processes  (C.(t)},  with  the 
processes  independent  for  different  values  of  i  and  having  uniform 
distributions  and  short  correlation  times,  then  the  phase  terms 
r^  sin  (w^  t  +  cp.)  would  have  no  effect  and  the  rotation  of  the  body 
would  be  immaterial.) 

Three  questions  of  interest  can  be  raised: 

a.  Suppose  the  radar  is  coherent  (observes  both  amplitude  and 
phase),  and  in  fact  forms  Bayes  estimates  predicated  on  the  assumption 
that  (C(t)}  =  0,  but  (C  (t)  }  is  in  fact  not  zero.  What  estimation 
errors  are  produced  by  {^(t)}? 

The  analytical  procedure  to  be  employed  in  answering  this  question 

is  as  follows.  The  estimates  i  or  0  are  the  same  as  those  given  in 

Section  II,  specialized  to  rotating  targets  as  in  Section  III-B, 

with  (C(t)}  assumed  to  be  zero.  Thus,  estimates  of  linear  forms  in 

a  are  still  linear  forms  in  S(t),  now  interpreted  to  include  additive 

noise  (s(t)},  and  estimates  of  quadratic  forms  in  a  are  still  quad- 

*»\ 

ratic  forms  in  S(t).  In  determining  the  moments  of  the  estimation 
error,  the  same  procedure  is  followed  as  in  Section  II. D,  except 
that  Eq.  (50)  can  no  longer  be  used.  In  place  of  this,  one  must  use 

S  =  F  {£}  CY  +  e  (50a) 

*+  w*  mm  U 

where  F  [C,]  is  formed  by  including  the  random  phase  process  Q  in  the 
factors  multiplying  each  component  of  a- 
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If  £(t)  is  bounded  (with  high  probability)  to  less  than  about 
one  radian  over  the  whole  observation  time,  then  F (£)  can  be  ex¬ 
panded  in  a  first  order  expansion 

F(r.,  cpi,  p ,  C,  t)  ^ 

F(r.,  cp.,  p,  0,  t)  +  C(0  Fc(r.,  cp.,  p,  0,  t)  (122) 

where  Fr  represents  a  derivative  in  an  obvious  way.  The  estimation 

t3 

errors  are  then  linear  or  quadratic  functions  of  both  (e(t))  and 
tC  (t) } ;  their  moments  can  be  determined  by  the  same  techniques  as 
in  Section  II. D.  In  such  a  case,  it  could  not  be  said  that  coherent 
processing  has  been  destroyed  by  {£(t)},  but  only  that  the  results 
have  been  degraded  by  an  additional  noise  term. 

For  coherence  to  be  destroyed,  it  would  be  necessary  to  assume 
that  [C(0]  is  not  bounded  to  a  small  value  over  the  observational 
interval.  If  [C(t)j  is  not  bounded  in  this  manner,  then  the  first- 
order  expansion  is  not  valid.  However,  it  is  still  possible  to  get 
the  first  and  second  moments  of  linear  or  quadratic  forms  in  S  in 
certain  cases,  e.g.,  if  (C(t)]  is  a  Gaussian  process  (defined  over 
-co,  co) ,  or  if  C(t)  has  a  correlation  time  which  is  short  compared 
with  the  time  in  which  r^  sin  (u)^_  t  +  cp^)  changes  significantly,  for 
all  i.  Methods  for  doing  this  form  another  entire  subject  and  will 
not  be  further  described  here. 

Returning  now  to  the  other  questions  of  interest: 
b.  Suppose  the  Bayes  estimate  is  formed,  utilizing  the  a  priori 
representation  Eq .  (121)  for  the  signal,  with  some  definite  a  priori 
distribution  attributed  to  {£(t)j.  What  are  the  estimates  then? 
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In  principle,  iC(t)]  can  be  regarded  as  a  parameter  vector  \  of  the 
type  treated  in  Section  II. B,  but  the  problem  still  remains  of 
elucidating  the  effect  of  this  particular  type  of  unknown  parameter 
vector.  This  in  turn  is  related  to  the  third  question: 
c.  Suppose  the  radar  only  observes  amplitude? 

The  Bayes  estimate  formed  by  taking  into  account  the  presence 
of  {£  (t) }  takes  the  following  form.  Let 


¥  =  FT  T)  F  =  B  -  A-1 


(123) 


Then,  from  Eq .  (121),  if  a'  is  large 


nr.,  <P.,  1,  r.,  cp.,  2) 


=  Y(r.,  cp.,  2,  r  ,  cp  ,  1)  =  0 


(124) 


and 


nr.,  cp.,  1,  r_.,  cp^,  1) 


=  Y(r. 

,  2, 

rr 

V 

2) 

•— 1 1  2° 

1! 

cos  1  r . 

L  1 

sin 

(”r 

t  +  cp.) 

-  r . 
J 

sin  (u! 

r 

t  + 

‘'j5’ 

(125) 

This  is  exactly  the  same  as  if  the  process  [C(b)}  were  identically 

T 

zero.  Moreover,  it  can  be  similarly  verified  that  F  A  F  is  inde¬ 
pendent  of  [C  (t) } ,  under  the  assumption  that  A  is  diagonal.  Thus, 

*  a  s 

the  matrices  B  or  q  entering  into  the  expression  for  a  or  p  are 
independent  of  (£(t)};  the  only  dependence  of  a  or  $  on  (C(t)} 
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arises  from  the  dependence  of  F  on  [£(t)},  in  the  expressions 
Eq.  (9)  or  (37).  To  be  specific,  let 


,  (1)  (2), 

o’  -  (a  ,  a  ) 


where  the  components  of  cf  ^  are  o^r.,  cp .  ,  1)  and  similatlv  for  a  ^  ^ 
Also,  B  consists  of  two  identical  blocks  corresponding  to  u  =  1  and 
p  =  2 ;  let  each  of  these  blocks  be  denoted  simply  B.  Also,  let 


B  (r  .  ,  cp .  ,  r  .  ,  cp . )  =  B.  . 
i  i  j  J  ij 


Then 


»"’«>  -Hr1).  -f. 

i  “  • 


>in  jjju^  t  +  r^  sin  (u^  t  +  +  C(t)  |  *  S(t) 


i<2)<o 

1  w 


ij 


2_ 

N 


cos  I  w  t  +  r.  sin  (to  t  +  o.)  +  C(t) 
L  c  l  r  l 


S(t) 


Now  further  assume  that  the  radar  is  a  pulse  radar,  subject  to 
the  following  conditions:  (1)  the  pulses  are  sinusoids  at  frequency 
and  have  length  much  greater  than  /.he  a  priori  range  extent  of  the 
target;  and  (2)  the  pulse  length  is  such  that  the  radial  motion 
of  all  a  priori  significant  scatterer  locations,  and  the  change  in 
£(t),  are  negligible  compared  to  a  wavelength  during  the  time  extent 
of  a  pulse.  (Radial  motion  means  motion  along  the  line  of  sight  to 


(126) 


(127) 


dt  (128 


dt 


the  radar ./ 
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The  matrix  B  can  be  found  from  Eqs .  (125)  and  (123);  in  this 
case,  Eq .  (125)  takes  the  form 


Y( 


r  i ,  cp^  1,  r.,  cp  ,  1)  =  ^(r^  cpi>  2,  r  ,  ,  2) 


=  )  cos  r.  sin  (id  t  +  cn.)  -  r.  sin  (id  t  +  cp.) 

N  L>  Li  rn'i  j  rn  l 

o  J 

n 

If  the  further  assumption  is  made  that  the  motion  (projected 
on  the  line  of  sight)  of  all  a  priori  significant  scatterer  locations 
is  negligible  compared  to  a  wavelength  during  the  interpulse  period 
At,  then  Eq .  (125a)  can  also  be  expressed 


Y(r.,  cp.,  1,  r.,  cp.,  1)  =  Y(r.,  cp.,  2,  r  ,  cp  2) 


1 

I  cos  |  r.  sin  (id  t  -f  cp.)  -  r.  sin  (u j  t  +  cp . ) 

J  L  l  r  i  j  r  j  _ 


N  At 
o 


dt 


which  is  the  same  as  Eq .  (125)  except  for  the  constant  factor  T  / At. 

P 

✓s 

To  see  what  form  a  takes,  let 


S(t)  =  S.(t)  sin  cju  t  +  S„(t)  cos  id  t 
1  c  2  c 


Then,  Eq .  (128)  becomes 


'y  _ 

S(1)(C)  =y(B1i..  cos  r.  sin  (U)  t  +  cp  , )  +  £ 

i  L>  /il  N  l  In  L  i  r  n  i  n 


o  n 


+  S„  sin  i  r.  sin  (ui  t  +  cp.)  +  r  i  >• 
2n  L  l  r  n  i  ^nJJ 


ij  h 


o  „  COS 
2n 


r,  sin  (w  t  +  cp . )  +  r  j 
.  l  r  n  l  nj 


S ,  i  ,  sin  'i  t  +  cp . )  C 

In  _  l  r  n  *n 


(125a) 


(125b) 


(129) 


(128a) 


3 
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where 


sm  '  W 


S2u  -  W 


C  -  Mt  ) 

n  n 


2  _  2 

Next,  consider  estimates  of  J  +  |  .  Let  be  the 

quadratic  form  associated  with  this  functional  of  a-  The  term 
Trace  ^  in  Eq .  (16)  is  independent  of  £.  The  second  term 

in  Eq .  (16)  becomes 


iT(C)  Q(l)  5(C) 


--{■i)  l  l  (t1)..  (?'*)  Vls2« 


cos  p .  -  S.  sin  p . 

jm  lm  jm  \ 


j  >k  m,n 


x  |  S„  cos  p.  -  S,  sin  p 


2n 


kn  In 


kn 


>)  j  i  n.  »■),{!. 

J ,k  m,n  J  k 


,  cos  p .  +  S„  sin  p  . 

lm  jm  2m  jm.J 


x  }_^in  cos  Pi.„  T  S0_  sin  p 


kn  2n 


knj 


where 


P.  =  r .  s i n  (dJ  t  +cp.)+C 
in  l  r  n  ri  sn 


This  simplifies  to 


(130) 


(131) 


(132) 
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ST(C)  Q(l)  a(C) 


T 

=  f  -E 

N 

o, 


I  I  vt1)  (l'1) 

i  j  ik 


j ,k  m, n 


L“lm  “In  '  “2m  “2nJ  LKjm 


r( 


S1„  ST„  +  S0„  S0_  I  COS  |_P  j 

Lpi= 


kn 


+  LSlm  S2n  +  S2m  Sl„ 


sm 


-  P, 


jm  kn 


-T 


Similarly,  the  expression  a  B  a  appearing  in  Eq .  (29a),  with 


X  =  b,  is 


(133) 


a  (O  B  a  (O 

2 


'  T 

=  I  _E 

N 

o 


I  Z  (t1)  {[Sl.  Sin  +  S2n,S 


2nj  cos  L^im  ^jn 


i,j  m , n 


+  _Slm  S2n  +  S2m  Sln 


sin  Ip.  -  P . 
L  im  i 


Jn. 


(134) 


Taking  account  of  Eq,  (132),  Eqs.  (133)  and  (134)  can  be  written 


Q' 


T(C)  o(C)  =  V*S.^  c.os  (c  -  C..  +  p 


L_j  ^mn 
m,  n 


(i)) 

m  n  mn  ' 


(133a) 


a1  (Q)  B  S(C)  =)  n 


cos 


m,n 


(134a) 


The  a  priori  p.d.f.  q (\ )  in  Eq .  (29a),  with  X  =  £,  can  be  assumed 
to  be  such  that  the  components  Q  of  C  are  mutually  independent  and 
uniformly  distributed  over  (0,  2rr)  •  The  evaluation  of  integrals  of 
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the  type  appearing  in  Eq .  (2.9a),  resulting  from  substitution  of 

Eqs.  (133a)  and  (134a)  into  Eq.  (29a),  seems  to  be  very  difficult  in 

general.  The  problem  clearly  amounts  to  evaluation  of  the  expected 

values  of  cos  cos  £  ,  sin  C  sin  C  ,  and  sin  C  cos  C  ,  with 
m  n  m  n  m  n 

respect  to  the  a  posteriori  distribution  of  j  implicitly  defined 
by  inserting  Eq .  (134a)  into  Eq .  (29a).  However,  there  are  some 
special  cases  in  which  such  evaluation  may  be  considerably  simplified. 
For  example,  suppose  the  a  posteriori  p.d.f.  of  Q  is  also  such 


rr 

that  jC 
t  n 

clear  that  this  will  not  in  general  be  true;  however,  it  will  tend  to 
be  true,  if  the  sampling  interval  At  is  such  that  a  significant  portion 
of  a  priori  scatterer  positions  move  along  the  line  of  sight  by,  say, 
a  wavelength  in  time  At. 


|  are  mutually  independent  and  uniformly  distributed.  It  is 


Then,  the  unconditional  estimate 
would  simply  be,  from  Eq .  (133), 


X  cos 


r .  sin  ('JU  t  +  cp.) 
I  J  r  n  Tj 


r.  sin  (U)  t  +  C+J.  )  | 
K  r  n  k  J 


(135) 


This  depends  only  on  the  square  amplitudes  of  the  observed 

signals . 

Equation  (135)  represents  the  estimate  achieved  when  fL(£)  is 
integrated  with  respect  to  the  a  priori  p.d.f.  of  It  is 

tempting  to  conjecture  that  this  represents,  in  general,  the  optimum 
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Bayes  estimate  of  0^  when  incoherence  is  due  to  the  fact  that  the 
radar  observes  only  the  signal  amplitudes.  When  the  components 
are  a  posteriori,  as  well  as  a  priori,  independent  and  uniformly 
distributed,  this  conjecture  must  be  true.  The  author  has  so  far 
been  unable  to  verify  the  truth  of  this  conjecture  in  general,  and 
in  fact  it  is  somewhat  doubtful  that  it  would  prove  to  be  true  in 
general.  It  is  easy  to  see  that  the  Bayes  estimate  when  the  radar 
observes  only  amplitude  must  always  be  at  least  as  good  as  the 
estimate  given  by  Eq .  (135). 


If  are  not  a  posteriori  independent  and  uniformly  dis¬ 

tributed  (as  will  be  the  case  when  the  sampling  rate  is  so  high  that 
a  significant  part  of  the  a  priori  scatterer  positions  move  along  the 
line  of  sight  by  less  than,  say,  one-quarter  wavelength),  then 
Eq .  (135)  does  not  represent  the  Bayes  estimate  if  the  radar  observes 
both  amplitude  and  phase.  There  are  still  cases  in  whicn  it  is 
probably  justified  to  assume  certain  simplifications.  For  example, 
suppose  the  sampling  rate  is  such  that  a  significant  part 
of  the  a  priori  scatterer  positions  move  along  the  line  of  sight  by 
the  order  of  a  wavelength  in  time  K  At,  K  >  1.  Then,  it  is  very 
probably  true  that  £  is  statistically  independent  cf  if 
|  m  -  n  |  >  K,  and  that  the  marginal  a  posteriori  p.d.f.  of  any  single 

Q  is  uniform  over  (0,2rr).  This  would  permit  one  to  set  equal  to 
n 

zero  all  terms  for  which  j  m  -  n|  =  K  in  integrating  Eq -  (133)  with 
respect  to  the  a  posteriori  p.d.f.  of  £. 


1 
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IV.  INITIAL  STEPS  TO  IMPLEMENTATION  OF  COMPUTER  PROGRAMS 


It  would  be  of  interest  to  implement  some  of  the  techniques  by 
programming  them  on  a  computer.  Two  approaches  are  indicated:  (a) 

A 

implementation  of  the  methods  for  forming  the  estimates  3,  and  (b) 
implementation  of  the  calculation  of  mean  square  error  for  linear 
and  quadratic  estimates. 

The  formation  of  the  estimates  themselves  would  be  a  Monte  Carlo 
program  since  it  would  be  necessary  to  produce  simulated  data  with 
additive  errors  or  phase  errors  (or  both).  Mean  square  error  could, 
of  course,  be  determined  by  repeated  trials  of  such  a  Monte  Carlo 
program.  Calculation  of  mean  square  error  via  the  formulas  of 
Section  II  would  not  require  Monte  Carlo  trials  or  generation  of 
simulated  data. 

It  would  be  most  practical  to  attempt  these  implementations  in 
an  ascending  order  of  complexity  with  respect  to  various  features  of 
the  problem,  such  as: 

Type  of  A  Priori  Knowledge  of  Scatterer  Distribution 

1.  Noncollinear ,  radially  symmetric  distributions  of  scatterer 

locations . 

2.  Collinear  distributions  without  a  specular  flash. 

3.  Collinear  distributions  with  a  specular  flash. 

4.  Mixture  of  (1)  and  (3),  or  (2)  and  (3). 


5. 


Surface  of  an  opaque  rotating  sphere. 
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Observation  Interval 

1.  Integral  number  of  rotation  periods,  via  Eqs .  (35)-(37),  (95), 

and  (100  ). 

2.  Fractional  period  observation  interval  sufficiently  short  no 

that  the  ratio  of  observation  time  to  the  time  required  for 
all  a  priori  significant  scatterer  locations  to  move  along 
the  line  of  sight  by  a  half-wavelength  is  a  small  integer. 

The  observation  times  could  then  be  considered,  with  good 
approximation ,  to  be  a  small  number  of  instants,  and  Eqs. 

* 

(35)- (37)  could  be  used  with  direct  matrix  inversion  to  find  h  • 

3.  Integral  number  of  rotation  periods,  with  total  number  of 

significant  a  priori  scattering  points  (separated  by  one-half 
wavelength)  a  small  number.  This  could  be  done  via  Eq .  (86) 

with  direct  matrix  inversion  to  find  B  ^  . 

4.  Fractional  observation  period  with  noncollinear  scatterer 

distribution,  via  Eqs.  (86)-(88)  to  compute  B  . 

Degree  of  Coherence 

1.  Completely  coherent. 

2.  Pulse  radar  with  complete  pulse-to-pulse  incoherence  due  to 

phase  errors  in  data  (but  with  radar  observation  of  received 
signal  phase),  via  Eqs.  (131) - (135)  and  (29). 

3.  Pulse  radar  with  partial  incoherence  due  to  correlated  phase 

errors  in  data. 

These  operations  can  be  done  for  a  variety  of  different  quadratic 
functionals  and  for  a  variety  of  assumptions  concerning  the  a  priori 
state  of  knowledge  and  the  true  state  of  affairs.  For  example,  it  is 
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of  interest  to  compare  the  accuracy  with  which  a  "true"  collinear 
scatterer  distribution,  and  the  aspect  angle  of  Eq .  (101),  can 
be  estimated,  both  when  the  a  priori  distribution  is  noncollinear 
and  when  the  a  priori  distribution  ic  collinear. 

It  is  of  particular  interest  to  attempt  the  evaluation  of  the 
integrals  appearing  in  Eq .  (29a),  when  Eqs.  (133a)  and  (134a)  are 


inserted.  This  will  yield  information  not  only  on  the  comparison 

r*  (l)~j2 

between  incoherent  and  coherent  processing  to  determine  a.  + 

-  2  ^  1 

I  3.  l  _ _ 1.1 _ _  -• _ _ 1  _1 _ _ 


,  but  also  on  the  accuracy  with  which  individual  phase  error 
differences  can  be  inferred,  even  if  the  phase  errors  are  independent 


from  pulse  to  pulse. 
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